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1.  INTRODUCTION 

Over  the  past  three  years  we  have  tried  to  inercasc  the  effectiveness  of  the  data  interpretation 
from  the  DMSP  SSIES  instrumentation.  We  have  pursued  two  major  paths  to  this  goal  One  is  the 
development  and  refincmen*  of  algorithms  for  producing  the  geophysical  parameters  and  another  is 
the  examination  of  data  ensuring  that  the  algorithms  capitalize  on  the  state  of  our  knowledge  about 
the  ionospheric  processes  being  described. 

The  following  report  describes  the  software  development  tasks  undertaken  to  make  the  data  re¬ 
duction  software  more  robust.  These  include  adaptations  to  changing  ion  composition  along  the 
spacecraft  track  and  to  changes  in  the  high  latitude  convection  pattern  resulting  from  magnetic  activ¬ 
ity.  We  have  also  adapted  algorithms  for  identifying  the  large  scale  features  in  the  ionospheric  convec¬ 
tion  pattern  from  a  single  satellite  pass  through  the  pattern.  I n  addition  to  the  algorithm  development 
and  refinement,  we  have  also  started  to  identify  reproducible  features  in  the  high  latitude  motion  that 
exist  when  the  interplanetary  magnetic  field  is  directed  northward.  We  expect  that  by  associating  mea¬ 
sured  quantities  with  specific  behavior  in  the  plasma  motion  we  will  he  able  to  develop  a  characteriza¬ 
tion  of  the  flow  that  is  comparable  in  quality  to  that  presently  available  during  southward  IMF.  The 
following  sections  document  the  work  performed  in  more  detail. 

2.  DMSPPOTMOD 

DMSP  Potential  Model  or  DMSPPOTMOD,  which  was  delivered  at  the  end  of  1989  with  support 
from  a  previous  contract,  was  the  first  major  analysis  program  used  for  classifying  the  high  latitude 
environment  based  on  SSIES  data.  Much  of  the  work  during  the  past  three  year  period  has  been  in 
exploiting  and  upgrading  the  analysis  developed  for  DMSPPOTMOD.  This  section  gives  an  overview 
of  the  scope  and  purpose  of  DMSPPOTMOD.  Further  details  arc  given  by  Ucclis  and  Hairston 
(1990). 

The  Defense  Meteorological  Satellite  Program  (DMSP)  is  a  series  of  polar  orbiting  satellites  de¬ 
signed  to  observe  the  weather  on  Earth  and  to  monitor  the  near-  Earth  space  environment  at  840  km. 
The  Special  Sensor — Ions,  Electrons  and  Scintillation  (SSIES)  package  which  has  flown  on  all  DMSP 
satellites  since  F8  was  launched  in  summer  1987,  provides  a  measure  of  the  velocity,  density,  tempera- 
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ture,  and  chemical  composition  of  the  thermal  plasma  in  the  upper  ionosphere.  One  of  the  instru¬ 
ments  on  SSIES  is  the  Ion  Drift  Meter  (IDM)  which  faces  into  the  direction  of  the  satellite  velocity 
and  measures  the  vertical  and  horizontal  components  of  the  ion  flow  perpendicular  to  the  satellite 
velocity.  Each  velocity  component  is  sampled  six  times  a  second.  Taking  the  perpendicular  flow  of 
the  ions  (v)  and  the  magnetic  field  at  the  satellite  (S),  the  electric  field  parallel  to  the  satellite’s  veloc¬ 
ity  vector  ( E )  can  be  calculated  using  the  equation 

E  =  -v  x  B. 

Integrating  the  electric  field  along  the  satellite  track  provides  a  measure  of  the  electrostatic  potential 
along  the  track. 

In  general,  the  flows  and  potentials  are  small  everywhere  along  the  satellite  orbit  except  in  the 
polar  regions.  Here  potential  differences  on  the  order  of  tens  of  kilovolts  exist  across  the  polar  region, 
fhe  magnitude  and  distribution  of  this  potential  are  important  parameters  for  understanding  the  ion¬ 
ospheric  environment.  Most  notably  here,  these  parameters  are  important  inputs  to  the  Magneto- 
spheric  Specification  and  Forecast  Model  (MSFM)  developed  by  Rice  University  for  use  by  the  Air 
Weather  Service.  DMSPPOTMOD  is  a  program  that  can  take  the  ion  flow  data  from  the  DMSP 
SSIES  package  and  produce  a  measurement  of  the  polar  cap  potential  that  can  be  used  by  the  MSFM. 

Data  for  one  half  orbit  going  from  an  equatorial  crossing  to  the  next  equatorial  crossing  are  averaged 
into  four  -  second  bins.  The  horizontal  flow  data  are  corrected  to  remove  the  effect  of  corotation.  The 
baselines  of  both  components  of  the  flow  arc  established  by  minimizing  the  flow  near  50  degrees  mag¬ 
netic  latitude.  The  50  degree  magnetic  latitude  is  assumed  to  be  at  zero  potential,  and  thus  under 
steady  state  condition  the  potential  should  be  zero  at  both  points  where  the  spacecraft  crosses  this  lati¬ 
tude.  However,  this  rarely  occurs  either  because  there  are  some  changes  in  the  potential  distribution 
during  the  twenty  minutes  it  takes  the  satellite  to  traverse  the  polar  region  or  because  there  exist  sys¬ 
tematic  offsets  in  the  measured  velocity.  For  most  polar  passes  an  offset  in  potential  remains  in  the 
calculated  potential  distribution  when  the  spacecraft  rccrosscs  the  50  degree  magnetic  latitude.  A  cor¬ 
rected  potential  distribution  which  has  both  ends  at  zero  potential  is  generated  by  taking  the  value  of 
this  offset  and  doing  a  linear  correction  to  the  measured  potential.  This  corrected  potential  distribu- 
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tion  is  used  to  determine  the  magnitude  and  locations  of  the  maximum  and  minimum  potential  seen 
along  this  pass,  as  well  as  the  location  of  the  zero  potential  point  between  the  maximum  and  minimum. 
The  geographic  asymmetry  of  the  location  of  the  zero  point  (i.c. — the  relative  distance  from  the  zero 
point  to  the  location  of  the  maximum  and  minimum)  is  used  to  designate  which  of  the  three  Heppncr- 
Maynard  patterns  (Hcppncr  and  Maynard,  1987)  best  fits  this  pass.  The  level  of  noise  in  the  flow  data 
inside  the  polar  cap  is  used  to  determine  a  quality  flag  for  the  data  for  this  pass  and  the  highest  magnet¬ 
ic  latitude  reached  by  the  satellite  is  used  to  determine  a  correction  factor  to  be  applied  to  the  ob¬ 
served  potential  maximum  and  minimum.  Finally,  these  parameters  along  with  a  unique  identifica¬ 
tion  number  for  the  pass  (SFINDEX)  are  saved  in  a  data  file  called  MSMDATA.DAT. 

3.  DMSPDBASE  AND  SUBSEQUENT  UPGRADES 

After  DMSPPOTMOD  was  delivered  to  Geophysics  Lab,  the  program  was  upgraded  in  spring 
1990  to  run  successive  passes  and  save  the  data  in  a  database  format.  This  new  version  was  renamed 
DMSPDBASE1  and  delivered  in  early  summer  1990.  Since  then  there  have  been  four  revisions  to 
the  program  with  the  latest,  DMSPD3ASE4,  being  delivered  in  January  1993  with  this  report.  Each 
revision  improved  on  the  quality  of  the  data  analysis  of  its  predecessor  and  was  based  on  our  increased 
understanding  of  the  data.  We  shall  examine  each  of  these  upgrades  and  its  results  in  turn. 

3.1  DMSPDBASE1 — Building  the  Database 

In  the  original  DMSPPOTMOD  program  the  values  for  the  ion  flow  velocity  averages,  their  stan¬ 
dard  deviations,  the  potentials,  and  the  satellite’s  locations  were  calculated  for  each  four -second 
timestep  during  a  polar  pass.  These  data  were  necessary  to  calculate  the  overall  potential  distribution 
during  a  pass  and  to  determine  the  parameters  that  were  saved  in  the  MSMDATA.DAT  data  file,  but 
once  these  parameters  were  found,  the  four-second  data  were  discarded.  With  increased  data  han¬ 
dling  capability  at  both  UT-Dallas  and  at  Phillips  Laboratory/GP,  this  four-second  data  could  also 
be  saved  for  the  database.  DMSPDBASE  1  adapted  the  DMSPPOTMOD  program  so  that  it  could 
read  an  entire  magnetic  tape  containing  ten  days  of  telemetry  data  in  the  Time -History -Database 
format.  As  it  processed  the  data  on  a  polar  pass  by  polar  pass  basis,  it  would  write  the  data  into  two 
files  which  are  referred  to  as  the  “shortfile”  (DMSPPARAM.DAT)  and  the  “longfile”  (DMSPPOT- 
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BIN. DAT).  The  shortfilc  contained  only  the  key  parameters  (potential  maximum  and  its  location, 
zero  crossing  location,  model  number,  etc.)  for  each  pass  and  was  identical  to  the  MSMDATA  file  gen¬ 
erated  by  DMSPPOTMOD.  The  longfile  contains  the  geophysical  data  about  the  IMF  and  Kp  indices 
before,  during,  and  after  the  pass  followed  by  the  measured  flow  data  and  calculated  potential  data 
for  each  four-second  period  covering  half  an  orbit  from  equatorial  crossing  to  equatorial  crossing. 
The  longfile  was  created  as  an  unformatted  sequential  File. 

The  format  for  the  longfile  is  given  in  the  following  tables.  The  first  table  gives  the  format  for 
the  first  line  of  each  longfile. 


PARAMETER 

VARIABLE  NAME 

DESCRIPTION 

index  number  for  pass 

SFINDEX 

unique  index  number  for  this 
polar  pass,  format  based  on 
satellite  number,  year,  day  of 
year,  hour,  and  minute  when 
this  pass  started  at  the  equator 
(NNYYDDDHHMM) 

IMF  Bx  component 

BXPREV 

these  arc  the  hourly  averaged 
IMF  components  during  the 
hour  PRIOR  to  start  of  pass 

IMF  By  component 

BYPREV 

IMF  Bz  component 

BZPREV 

IMF  Bx  component 

BXSTRT 

these  are  the  hourly  averaged 
IMF  components  during  the 
hour  AT  the  start  of  the  pass 

IMF  By  component 

BYSTRT 

IMF  B7  component 

BZSTRT 

IMF  Bx  component 

BXSTOP 

these  are  the  hourly  averaged 
IMF  components  during  the 

hour  AFTER  the  start  of  pass 
(Note  that  aii  the  IMF  values 
are  given  in  nanoteslas  in  the 
geocentric  solar  magnetic 
(GSM)  coordinate  system 

IMF  By  component 

BYSTOP 

IMF  Bz  component 

BZSTOP 

AE  index  at  start  of  pass 

IAEINDEX1 

not  yet  available,  fill  data 

AE  index  at  end  of  pass 

IAEINDEX2 

not  yet  available,  fill  data 

Kp  index  at  start  of  pass 

KPSTART 

Kp  index  is  coded  as  a  2 -digit 

integer  as  follows:  Kp  0  =  0, 
KpO-t-  =  1 ,  Kp  1  -  =2,  Kpl  = 

3,  Kp  1  +  =4, . Kp  9  -  =  26, 

Kp  9  =  27 

Kp  index  at  end  of  pass 

KPSTOP 

always  equal  to  810,  or  54  min¬ 
utes  of  data,  serves  as  a  check 
on  the  size  of  the  file 


number  of  4-second  data 
blocks 


1MAX 


horizontal  flow  correction 

CHF 

vertical  flow  correction 

CVF 

offset  applied  to  the  horizontal 

and  vertic  il  flow  data  in  order 
to  bring  the  endpoints  to  zero 
velocity 


(Note  that  for  a  given  pass,  theSFINDEX  number  of  the  longfilc  is  the  same  as  theSFINDF.X  number 
of  the  shortfile  for  that  pass.  This  allows  for  easy  correlation  of  passes  between  files  in  the  two  data¬ 
bases.  For  the  hourly  averaged  By,  By,  and  B,  components  of  the  IMF.  note  that  in  this  analysis,  the 
hour  in  which  the  pass  begins  is  defined  as  the  hour  when  the  satellite  crossed  the  equator,  even  if  that 
occurs  during  the  59th  minute  of  that  hour.  If  there  arc  no  IMF  data  available  for  any  of  these  times. 


then  fill  values  of  999  arc  saved  in  their  place  to  the  File.  The  rationale  for  recording  this  IMF  data 
is  to  allow  a  search  of  the  database  for  passes  which  occurred  during  times  when  the  IMF  was  steady 
over  a  period  of  two  or  three  hours.  As  the  AE  indices  are  not  yet  available,  Fill  data  of  990  is  saved 
to  the  file.  At  some  point  in  the  future  when  this  data  becomes  available,  it  will  be  incorporated  into 
the  database.  The  number  of  four-second  blocks  in  this  pass  is  given  as  the  value  of  IMAX.  This  is 
always  equal  to  810,  corresponding  to  54  minutes  worth  of  data.  Originally  the  number  of  four-sc- 
cond  blocks  for  each  pass  was  going  to  vary  depending  on  exactly  the  time  between  the  equatorial 


crossings,  thus  causing  IMAX  to  vary.  It  was  decided  to  make  all  the  passes  a  uniform  length  to  simpli¬ 
fy  the  processing  arid  the  IMAX  value  was  retained  as  a  check  on  this.  A  length  of  54  minutes  takes 
each  of  the  satellites  from  an  equatorial  crossing  to  slightly  past  the  next  equatorial  crossing.  This 
results  in  some  overlap  of  data  in  the  files,  but  this  is  preferable  to  the  possibility  of  missing  some  data 
by  choosing  too  short  a  time  length.  The  final  values  in  the  arc  the  values  for  the  horizontal  flow 
correction  and  the  zerttcal  flow  correction  which  were  applied  to  the  de -corotated  flow  data  to  bring 
the  endpoints  of  the  flows  for  this  pass  to  zero  (CHF,  CVF).  These  are  saved  for  later  work  in  checking 
for  systematic  variations  in  the  drift  of  the  baseline  of  the  flow  data.) 


The  next  tabic  gives  the  format  for  the  remaining  810  lines  of  data  in  a  longFile: 


PARAMETER 

VARIABLE 

time 

XUTIME 

horizontal  ion  flow  velocity 

FLWH3 

vertical  ion  flow  velocity 

FLWV3 

DESCRIPTION 


given  in  seconds  ( 0-86400 ) 


four  second  averages  given  in 
kilomcters/sccond  (plus  3 
km/s,  see  note  below) 
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standard  deviation  of  horizon¬ 
tal  ion  flow  velocity 

STDEVH 

standard  deviation  of  the  data 
points  included  in  the  four  se¬ 
cond  averages 

standard  deviation  of  vertical 
ion  flow  velocity 

STDF.VV 

number  of  points 

\MBPTS 

number  of  points  included  in 
the  four  second  averages,  nom¬ 
inally  equal  to  24,  the  number 
is  the  same  for  both  horizontal 
and  vertical  analysis 

clectrosta  potential 

POTLNG 

corrected  potential  given  in  ki¬ 
lovolts 

magnetic  latitude  {at  satellite) 

SCCHMIjVTS 

given  in  degrees 

magnetic  latitude  (at  110  km) 

SCCHMLAT 

given  in  degrees 

magnetic  local  time  (at  110 
km) 

SCCHMLTM 

given  in  hours  and  decimal 
fractions  of  an  hour 

invariant  latitude 

SCCH1NVLAT 

given  in  degrees 

geographic  latitude 

SCCHLAT 

given  in  degrees 

geographic  east  longitude 

SCCHLONG 

given  in  degrees 

(The  values  for  the  ion  flow  velocities  have  had  the  corotation  of  the  ionosphere  and  the  flow  offset 


values  (CHF,  CVF)  removed,  and  thus  are  given  in  the  inertial  frame  of  reference  fixed  to  the  sun 
-Earth  line.  In  order  to  consent-  the  space  a  negative  sign  would  occupy  in  the  record,  all  the  values 
for  the  flow  data  have  3.0  km/s  added  to  them  before  they  are  saved  to  the  data  file.  This  insures  that 
all  the  flow  values  arc  positive.  This  excess  is  always  removed  when  the  data  is  read  out.  The  standard 
deviation  and  number  of  points  in  the  analysis  arc  saved  as  a  check  on  the  noisiness  and  completeness 
of  the  data  incorporated  into  the  four-second  flow  averages.  The  values  for  the  electrostatic  poten¬ 
tial  are  'hose  of  the  corrected  potential  distribution  described  above.  The  geophysical  coordinates 
every  four  seconds  arc  interpolated  using  the  one  minute  increments  taken  from  the  orbit  data  in  the 
telemetry  ) 


The  blocks  of  data  for  each  pas^.  are  written  sequentially  onto  the  file  until  the  full  ten  days  of  data 
have  been  analyzed.  The  writing  of  this  data  to  the  longfilc  occurs  in  BLOCK  4B  of  DMSPDBASF1 , 
However,  if  the  pass  is  determined  to  be  unusable  for  some  reason,  then  a  null  file  of  zeroes  is  written 
to  the  shortfile  under  the  SFINDF.X  number  for  that  pass,  but  no  corresponding  longfilc  is  created. 
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3.2  DMSPDBASE1 B — Upgrading  the  Endpoint  Selection  Routine 

The  first  revision  came  later  in  summer  of  1990.  Asa  result  of  examining  the  data  from  the  data¬ 
base,  it  became  apparent  that  the  algorithm  for  choosing  the  location  of  the  endpoints  of  the  pass 
should  be  improved.  The  original  algorithm  placed  the  endpoints  symmetrically  about  the  center  of 
the  magnetic  coordinate  system.  Thus  both  ends  would  be  at  the  same  magnetic  latitude.  In  the  origi¬ 
nal  algorithm  both  ends  would  be  checked  at  50  degrees  magnetic  latitude.  If  either  point  was  unsuit¬ 
able,  then  the  algorithm  checked  both  ends  at  45  degrees  magnetic  latitude.  This  process  repeated 
until  suitable  endpoints  were  found  or  the  search  terminated  at  20  degrees  magnetic  latitude.  While 
this  worked  for  the  majority  of  the  cases,  sometimes  passes  were  encountered  in  which,  for  example, 
the  starting  endpoint  at  50  degree  magnetic  latitude  would  he  acceptable  while  the  stopping  endpoint 
was  not.  The  program  would  check  at  45  degrees  magnetic  latitude  and  then  find  the  new  starting  point 
was  unsuitable,  but  the  new  stopping  endpoint  was  acceptable.  In  such  a  case  the  obvious  solution 
would  be  to  use  the  starling  endpoint  at  50  degrees  and  the  stopping  endpoint  at  45  degrees.  Instead 
the  original  algorithm  would  move  on  to  the  40  degree  endpoints  and  check  those  points.  At  best, 
both  those  points  would  be  chosen,  giving  the  pass  an  excess  “tail”  on  either  end.  or,  as  happened  in 
some  cases,  no  case  was  found  where  both  endpoints  at  the  same  latitude  were  suitable.  This  resulted 
in  a  usable  pass  being  ruled  as  unsuitable  and  discarded  from  the  database . 

The  modification  was  fairly  simple.  Block  2A  was  rewritten  such  that  the  search  routine  consid¬ 
ered  each  endpoint  separately.  If  the  starting  endpoint  at  50  degrees  was  suitable  but  the  stopping 
endpoint  was  not,  then  the  program  only  continued  searching  for  a  suitable  stopping  endpoint.  Once 
that  was  found,  then  the  rest  of  the  analysis  on  the  pass  continued. 

3.3  DMSPDBASE2 — Modifying  the  Highest  Latitude  Passes 

As  more  passes  were  processed  it  was  realized  that  errors  in  determining  the  ground  track  of  the 
satellite  in  a  magnetic  coordinate  system  were  occurring  frequently  during  the  passes  that  went  above 
80  degrees  magnetic  latitude.  The  two  problems  addressed  here  were  an  error  in  interpolating  the 
four-second  increments  of  the  spacecraft’s  location  in  magnetic  coordinates  and  an  error  in  deter- 
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mining  the  magnetic  latitude  above  85  degrees.  noth  these  problems  were  fixed  in  spring  1991  and 
the  updated  version  after  this  was  renamed  DMSPDBASH2. 

In  the  original  DMSPDBASEI  program  the  four -second  increments  of  magnetic  latitude,  mag¬ 
netic  local  time,  invariant  latitude,  geographic  latitude  and  geographic  longitude  were  determined  by 
a  simple  linear  interpolation  of  the  values  given  every  minute  in  the  orbit  data.  This  worked  satisfacto¬ 
rily  in  most  cases,  but  failed  at  latitudes  higher  than  80  degrees.  The  failure  arises  from  treating  polar 
coordinates  as  if  they  were  rectilinear  coordinates  in  the  interpolation  and  simply  interpolating  each 
component  individually.  At  lower  latitudes  where  there  was  not  much  change  in  the  magnetic  local 
time  in  the  one -minute  intervals,  thus  very  little  error  resulted  from  using  such  a  method.  But  at 
the  highest  latitudes  where  there  was  a  fairly  large  change  in  the  magnetic  local  time  over  the  course 
of  one  minute,  this  caused  “scallop  -shaped”  curves  to  appear  in  the  satellite's  track  (Figure  1 ).  This 
was  corrected  by  changing  the  interpolation  procedure  in  the  COROTFIX  subroutine  where  the 
four-second  increments  were  calculated  and  renaming  the  subroutine  COROTF1XG.  In  the  new 
interpolation  procedure,  the  polar  coordinates  arc  transformed  into  rectilinear  coordinates  and  then 
a  linear  interpolation  is  performed  on  them  in  order  to  obtain  the  four- second  increments.  After  this, 
the  new  rectilinear  coordinates  of  the  four-second  increments  arc  transformed  back  into  the  original 
polar  coordinate  system,  thus  giving  a  straight  ground  track  in  the  polar  coordinate  system. 

The  second  problem  occurred  only  when  the  satellite  went  above  85  degrees  magnetic  latitude. 
Magnetic  latitude  is  computed  by  tracing  the  field  line  passing  through  the  satellite  down  to  1 10  km 
and  measuring  the  angle  between  that  point  and  the  magnetic  dipole.  The  model  used  to  determine 
the  magnetic  latitude  for  the  orbit  data  in  the  DMSP  header  tended  to  have  large  errors  in  the  calcula¬ 
tions  the  closer  the  magnetic  latitude  got  to  90  degrees.  Thus  on  some  passes  that  got  close  to  the 
magnetic  pole,  the  one -minute  locations  in  the  magnetic  coordinate  system  showed  an  erroneously 
jagged  path  (Figure  2a)  where  one  to  three  locations  arc  pushed  away  from  the  true  ground  track. 
This  problem  was  solved  by  adding  Block  ID  to  the  program  which  takes  any  pass  that  goes  above 
85  degrees  magnetic  latitude  and  relocates  the  point(s)  above  85  degrees  using  a  third  order  polyno¬ 
mial  least  squares  fit  to  the  points  surrounding  it.  While  not  all  passes  which  went  above  85  degrees 


Magnetic  Local  Time  and  Magnetic  Latitude 


Figure  i .  This  plot  shows  the  ground  track  ot  a  DMSP-  FX  polar  pass  in  magnetic  latitude  and  mac 
netic  local  time  coordinates  using  the  incorrect  path  interpolation  from  the  original  DMSPPOT 
MOD  routine.  Note  the  small  “scallop”  shaped  humps  near  the  poles  that  result  from  this  interpola 
tion.  The  cusps  between  the  ’  scallops  arc  the  locations  ot  the  one  minute  marks  in  the  telemetry 
header.  Note  also  that  at  lower  latitudes,  this  interpolation  works  quite  well. 


0 

Magnetic  Local  Time  and  Magnetic  Latitude 


Figure  2a.  This  polar  dial  show  s  an  example  of  a  high  latitude  pass  with  a  single  point  out  of  align 
ment.  The  asterisks  denote  the  locations  of  the  satellite  at  the  one  minute  marks  given  in  the  teleme 
try  header. 


had  the  high  latitude  point  mislocatcd,  the  error  occurred  frequently  enough  that  it  was  deemed 
worthwhile  toappiy  the  routine  automatically  to  all  of  them.  The  rationale  for  this  was  that  if  the  high¬ 
est  latitude  point  was  not  mislocatcd,  then  the  smoothing  routine  had  very  little  effect  on  it 

The  smoothing  was  done  by  taking  the  locations  in  magnetic  latitude  and  magnetic  local  time  of 
the  seven  points  at  the  highest  magnetic  latitude  and  converting  them  into  a  rectilinear  coordinate 
system.  The  seven*  components  and  the  seven  y  components  along  with  their  corresponding  times 
arc  independently  run  through  a  least  squares  subroutine  (LESSQ3)  which  returns  a  set  of  coefficients 
for  a  third-order  polynomial.  The  program  calculates  a  new  set  of*  andy  coordinates  for  the  seven 
locations.  This  gives  a  corrected  location  for  the  erroneous  point  above  85  degrees  with  only  a  mini¬ 
mum  adjustment  to  the  locations  of  the  surrounding  points  (Figure  2b).  The  locations  for  the  new 
points  are  then  transformed  back  into  magnetic  latitude/magnetic  local  time  coordinates  and  replace 
the  original  values  for  of  magnetic  latitude  and  magnetic  local  time  in  the  arrays  SCMLAT  and 
SCMLT.  Note  that  this  correction  is  done  before  the  program  calls  the  COROTFIXG  subroutine  and 
thus  the  interpolated  four-second  increments  of  magnetic  latitude  and  magnetic  local  time  are  based 
on  these  smoothed  values,  not  the  original  erroneous  ones. 

3.4  DMSPDBASE3— Correcting  for  the  “H+  fuzz” 

Initially  the  algorithm  we  designed  to  perform  the  potential  analysis  assumed  that  the  vertical  and 

horizontal  ion  drift  were  close  to  zero  in  the  region  equatorward  of  the  auroral  zone.  (In  most  cases 
the  50  degree  magnetic  latitude  points  arc  chosen  ac  the  starting  and  stopping  point  and  the  flow  data 
are  close  to  zero  there.  The  algorithm  does  allow  for  other  points  to  be  chosen  if  the  polar  cap  region 
has  expanded,  as  in  the  case  of  the  “Great  Storm”  of  March  1989.)  The  baseline  for  the  entire  pass 
is  then  zeroed  using  the  average  value  of  these  endpoints  and  the  analysis  is  performed  on  the  cor¬ 
rected  flow  data.  During  the  period  of  December  1987  through  March  1988  an  anomaly  appeared  in 
the  data  of  DMSP-F8;  both  the  horizontal  and  vertical  flow  data  showed  large  values  (>  0.9  km/s) 
and  a  large  scatter  ( —0.4  km/s)  during  the  northbound  leg  of  the  orbit  at  MLT  06(H).  These  “humps" 
in  the  data  began  at  midlatitudcs  in  the  southern  hemisphere  and  continued  until  just  before  the  satel¬ 
lite  reached  the  northern  aurora!  region  (Figure  3).  The  “humps"  were  named  “H  +  fuzz”  and  analy- 
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Magnetic  Local  Time  end  Magnetic  Latitude 

Hi-'uie  2b.  Tii i s  polar  dial  shows  the  Mine  polar  pass  shown  in  figure  2a  attei  the  smoothing  al 
nthm  has  been  applied  to  the  misaligned  point  and  the  six  points  surrounding  it. 
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Figure  3.  This  plot  shows  the  number  density  and  plasma  I  low  data  tor  a  single  orbit  ot  I  S  on  1 
March  1988.  The  plots  starts  as  the  satellite  crosses  the  equator  heading  northward,  goes  cnei  the 
northern  polar  region,  past  the  equator  at  the  midpoint  ot  the  plot,  over  the  southern  polar  rest  ion. 
and  ends  at  the  next  northbound  equator  crossing.  The  H+  hi/./'  is  apparent  as  the  non  -  zero  velocities 
in  both  the  horizontal  and  vertical  I  low  s  on  either  side  of  the  northern  polar  region.  Nomin.illv  the 
How  velocities  in  those  regions  should  he  close  to  zero  such  as  are  seen  surrounding  the  southern 
polar  regions. 


sis  showed  that  they  were  an  artifact  of  the  IDM  whenever  the  light  ions  (hydrogen  and  helium  ions) 
concentration  was  greater  than  20  or  25%  of  the  total  ion  density.  This  is  consistent  with  the  condi¬ 
tions  seen  in  the  polar  ionosphere  during  that  hemisphere’s  winter,  particularly  on  the  nightside.  Since 
the  F8  orbit  is  tilted  at  98.8  degrees,  it  is  always  on  the  nightside  of  the  0600-1800  hours  local  time 
line  in  the  northern  hemisphere  and  on  the  daysidc  of  the  same  line  in  the  southern  hemisphere.  This 
explains  why  this  phenomenon  was  not  observed  in  the  southern  hemisphere  when  F8  was  first 
launched  in  June  1987  during  the  southern  winter.  Instead  it  First  appeared  only  in  the  northern  hemi¬ 
sphere  in  the  late  fall  of  1987.  The  H+  fuzz  did  not  appear  on  the  southbound  leg  of  the  northern 
hemisphere  orbit  because  that  occurred  in  the  early  evening  portion  of  the  ionosphere  at  MLT  1800. 
The  plasma  there  had  just  moved  in  from  the  daysidc  and  was  still  hot  enough  that  the  scale  height 
of  the  0  +  ions  was  large  and  the  relative  abundance  of  0+  at  800  km  was  >  75  -  80%>.  As  this  plasma 
moved  around  into  the  nightside,  it  cooled  off  and  the  scale  height  decreased.  By  the  time  the  plasma 
iiad  moved  around  to  0600  MLT  many  of  the  heavier  0+  ions  had  sunk  to  lower  altitudes  leaving  the 
lighter  H+  and  Hc+  ions  to  make  up  a  larger  proportion  of  the  plasma  at  800  km.  This  explains  why 
the  H+  fuzz  was  so  much  smaller  on  the  northern  duskside  which  had  just  recently  entered  darkness. 
In  the  southern  hemisphere  no  fuzz  was  seen  on  the  duskside  leg  of  the  orbit  because  that  portion  was 
still  sunlit.  On  the  southern  dawnside  the  H  4-  fuzz  would  appear  at  midlatitudcs  because  this  plasma 
had  only  recently  moved  into  the  sunlit  region. 

Curiously,  the  phenomenon  of  the  H  +  fuzz  did  not  reappear  in  the  next  winter  period  (December 
1988-March  1989)  nor  in  any  of  the  subsequent  winter  periods  through  March  1992  in  the  F8  data. 
We  believe  this  to  be  caused  by  the  increased  solar  activity  which  began  in  1988.  During  solar  maxi¬ 
mum  the  atmosphere  is  hotter  than  during  solar  minimum  and  in  turn  the  scaie  heights  of  the  various 
species  in  the  plasma  are  larger,  thus  keeping  the  proportion  of  0+  at  800  kilometers  above  the 
75  -80%  level.  While  wc  do  not  expect  this  H  +  anomaly  to  reappear  before  the  end  of  F8’s  operation¬ 
al  lifetime,  it  may  appear  in  the  FI  1  (and  subsequent  satellites  in  the  0400-  1600  MLT  orientation) 
data  in  the  next  few  years. 
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The  H+  fuzz  is  a  concern  because  the  original  algorithm  used  in  DMSPPOTMOD  and  DMSP- 
DBASE  tended  to  chose  starting  points  inside  the  H+  fuzz  region  for  passes  in  the  northern  hemi¬ 
sphere.  This  resulted  in  an  incorrect  baseline  being  set,  which  in  turn  produced  bad  results  for  the 
calculated  potential.  The  algorithm  in  DMSPDBASE  was  upgraded  in  the  fall  of  1991  so  that  it  would 
avoid  choosing  starting  points  inside  the  H  +  fuzz  regions.  Figures  4a  and  b  show  the  same  pass  before 
and  after  the  upgrade  of  this  algorithm  in  Block  2A. 

As  shown  in  these  figures,  many  passes  show  the  H  +  fuzz  extended  above  the  50  degrees  magnetic 
latitude  on  the  dawn  side.  For  this  reason  the  new  algorithm  starts  looking  for  a  suitable  starting  end¬ 
point  at  65  degrees  magnetic  latitude.  The  criteria  for  choosing  a  suitable  endpoint  was  determined 
empirically  by  examining  the  F8  flow  data  during  the  period  from  December  1987  through  March 
1988.  Anewsubroutinc  named  SEARCH  was  incorporated  into  the  program  which  used  the  following 
procedure  for  determining  whether  a  trial  point  at  a  given  latitude  is  a  suitable  endpoint: 

1)  It  checks  the  standard  deviation  for  the  four-second  horizontal  flow  average  (STDEVH)  and 
the  four-second  vertical  flow  average  (STDEVV)  as  well  as  the  standard  deviation  of  the  aver¬ 
ages  for  each  component  on  either  side  of  this  trial  point.  If  any  one  of  these  six  standard  devi¬ 
ations  is  greater  than  0.14,  then  the  data  at  this  location  is  deemed  to  be  too  noisy  to  be  suitable 
as  an  endpoint. 

2)  It  checks  how  far  the  averaged  flow  values  at  this  trial  point  deviate  from  zero.  For  most  of  the 
time  for  all  the  DMSP  satellites  the  horizontal  flow  should  be  within  ±0.2  km/s  of  zero  and  the 
vertical  flow  should  be  within  ±0.35  km/s  of  zero.  However  there  was  a  large  offset  in  the  flow 
baseline  in  the  early  months  of  F8,  so  for  the  1987  period  values  of  ±0.4  km/s  for  the  horizontal 
flow  and  ±  0.55  km/s  for  the  vertical  flow  are  used.  If  either  the  horizontal  flow  or  vertical  flow 
at  this  location  is  outside  these  limits,  then  this  point  is  deemed  unsuitable  as  an  endpoint.  This 
removes  the  chance  that  a  point  on  the  “hump”  of  the  H  +  fuzz  region  would  be  selected  as  an 
endpoint. 
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Kigure  4a.  This  plot  shows  the  analysis  of  the  electrostatic  potential  for  the  northern  polar  pass 
shown  in  figure  3  without  using  the  corrections  incorporated  into  PMSPDB ASF.3.  The  program 


erroneously  chose  the  45  degree  magnetic  latitude  point  on  the  daw  aside  w  Iir  it  is  inside  the  1 1+  \\i// 
region  (note  the  change  in  scales  for  the  flow  voles  ities  between  figures  3  and  4a).  The  difference 
between  the  flow  velocities  at  the  starting  and  stopping  points  are  about  1  kiloineter/second.  thus 


the  program  takes  the  average  ot  the  two  ends  and  sets  a  f  lawed  baseline  lor  the  analysis.  This  results 
in  the  potential  curve  ha\  ing  a  spurious  negative  potential  portion  between  45  and  hi 1  degrees  on  the 
daw  nside  which  w  recks  the  oseiall  shape  of  the  potential  curve. 
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Figure  4b.  This  plot  shows  the  same  pass  as  figure  4a  analyzed  using  the  new  algorithms  in 
DMSPDBASE3.  The  new  version  of  the  analysis  correctly  selects  usable  endpoints  at  62  and  60 


degrees  magnetic  latitude,  which  results  in  an  accurate  calculation  of  the  potential  curve. 


3)  It  takes  the  average  of  the  trial  point  along  with  the  two  four-second  averaged  flow  values 
prior  to  and  two  averages  following  the  trial  point.  It  calculates  the  variance  for  these  five  values 
and  checks  if  it  is  greater  than  0.006.  This  procedure  is  repeated  for  vertical  flow  values  around 
the  trial  point.  If  either  of  these  variances  is  greater  than  0.006,  then  the  region  around  this  point 
has  too  much  scatter  and  is  deemed  unsuitable  for  an  endpoint.  This  check  guards  against  choos¬ 
ing  an  endpoint  in  the  H+  fuzz  region  where  there  is  a  large  scatter  in  the  four-second  averaged 
flow  data. 

If  all  three  of  these  checks  are  passed,  then  the  SEARCH  subroutine  returns  a  value  of  1  for  the 
variable  IGO  to  indicate  that  this  point  is  suitable  as  an  endpoint.  The  selection  of  the  exact  values 
of  the  cutoffs  in  the  three  checks  in  SEARCH  was  performed  empirically  to  eliminate  the  chance  of 
accidentally  choosing  an  endpoint  in  the  H  +  fuzz  region.  These  values  may  be  changed  to  provide 
a  looser  or  tighter  selection  criteria. 

The  endpoint  algorithm  in  Block  2A  was  revised  to  change  the  order  of  the  search  for  endpoints. 
The  revised  routine  begins  by  checking  the  trial  starting  points  at  65,  60,  55,  and  50  degrees.  If  any 
of  these  were  suitable,  the  suitable  point  lowest  in  latitude  is  picked  as  the  starting  point.  If  none  of 
the  four  are  suitable,  the  algorithm  continues  to  inspect  trial  points  at  five  degree  increments  until 
either  a  suitable  starting  point  was  found  or  the  algorithm  reaches  20  degrees  latitude.  If  no  suitable 
starting  point  has  been  reached  by  20  degrees,  the  program  starts  the  search  again  at  64  degrees  and 
examines  trial  points  in  one  degree  increments.  This  continues  until  either  a  suitable  starting  point 
is  found  or  it  reaches  20  degrees  latitude.  If  the  algorithm  fails  to  find  a  suitable  starting  point,  the 
program  branches  to  where  it  writes  a  null  file  for  this  pass  to  the  shortfile  in  the  database,  then  cycles 
back  to  start  on  the  next  pass.  If  the  routine  docs  find  a  suitable  starting  point,  it  repeats  the  same 
procedure  to  find  a  suitable  stopping  point.  Doing  the  search  in  one  degree  increments  will  find  a  suit¬ 
able  endpoint  for  most  cases,  but  it  is  a  time-consuming  procedure  if  attempted  on  the  majority  of 
non  -  H  +  fuzz  cases.  Thus  the  rationale  behind  this  algorithm  is  to  try  to  five  degree  increment  search 
first  on  all  cases,  only  attempting  the  one  degree  increment  search  on  those  that  failed. 
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The  improvement  of  the  new  algorithm  in  DMSPDBASE3  is  shown  by  comparing  the  results  of 
the  earlier  version  to  the  revised  version  for  a  typical  period  when  the  H+  fuzz  appeared  For  the 
period  of  March  1  through  10, 1988  DMSP-F8  data  has  134  complete  northern  passes  and  1 33  south¬ 
ern  passes,  these  being  defined  as  periods  of  uninterrupted  data  from  equator  crossing  to  next  equator 
crossing.  There  is  a  data  gap  of  about  five  orbits  during  the  latter  half  of  day  65  through  the  beginning 
of  day  66.  Those  lost  passes  arc  not  included  in  the  totals  below.  For  this  comparison  each  polar  pass 
is  placed  in  one  of  four  categories: 

nominal — the  baseline  correction  and  choice  of  endpoints  was  done  properly,  no  "eyeball'' 
correction  would  improve  the  results  significantly.  Note  that  this  does  not  necessarily 
mean  that  the  data  are  good,  there  is  still  the  possibility  that  the  pattern  changed  during 
the  pass,  thus  giving  a  confused  or  garbled  result.  “Nominal”  means  that,  given  the  data 
for  this  particular  pass,  this  analysis  is  probably  as  good  as  it  will  ever  be. 

bad — cases  where  the  program  has  incorrectly  chosen  an  endpoint  in  either  the  “H  +  fuzz”  or 
the  auroral  region,  resulting  in  a  skewed  baseline  and  an  erroneous  potential  curve. 

data  gap — throughout  the  DMSP  data  stream  there  are  data  gaps  ranging  from  one  minute  to 
hours.  In  some  cases  where  the  gap  is  only  a  minute  or  two  in  length  or  else  occurs  out¬ 
side  of  the  polar  auroral  region  it  is  possible  to  recover  the  potential  curve  for  that  pass. 
However,  this  is  only  done  using  the  program  DMSPDBASEINDV  which  requires  a 
human  in  the  loop  to  decide  whether  the  pass  is  doable  or  not.  For  the  automatic  rou¬ 
tines,  neither  DMSPDBASE2  nor  DMSPDBASE3  has  the  sophistication  to  recover 
these  passes.  As  a  result  the  pass  with  the  data  gap  returns  a  null  file  to  the  database 
and  the  subsequent  pass  is  skipped  over.  Since  such  gaps  are  relatively  infrequent  (here 
<  2%)  this  is  not  considered  a  major  flaw  in  the  program. 

failure — a  case  where  the  new  algorithm  has  examined  the  data  from  65  to  20  degrees  magnetic 
latitude  on  both  sides  of  the  pass  using  steps  of  one  degree  and  was  unable  to  locate 


either  a  suitable  starting  or  stopping  endpoint.  Because  this  is  part  of  the  new  algo¬ 
rithm,  this  ease  does  not  exist  for  the  analysis  using  the  old  algorithm. 


Using  the  data  from  March  1  —  10, 1988  (day  61  -70)  the  old  and  new  algorithms  give  the  following 
results: 

Using  DMSPDBASE2 : 

North:  nominal  12(9%)  bad  120  (89.5%)  data  gap  2  (1.5%) 

South:  nominal  131  (98.5%)  bad  0(0%)  data  gap  2  (1.5%) 


Using  DMSPDBASE3: 

North:  nominal  128  (95.5%)  bad  1  (0.7%;) 
South:  nominal  131  (98.5%)  bad 0(0%) 


data  gap  2  ( 1 .5%)  failure  3  (2.3%) 
data  gap  2  ( 1 .5%)  failure  0  (0%) 


Obviously,  the  new  algorithm  vastly  improved  the  analysis  on  northern  passes  where  the  H  +  fuzz 
was  a  problem.  The  two  algorithms  gave  identical  results  in  the  southern  hemisphere  where  the  H  + 
fuzz  was  absent. 


3.5  DMSPDBASE4 — Upgrading  the  Heppner-Maynard  Model  Identification  and 
Further  Corrections  to  the  Endpoint  Selection  Routine 

From  the  summer  of  1992  through  December  1992  work  was  done  on  the  latest  version  of  the 
program  which  was  called  DMSPDBASE4.  Two  major  problems  were  addressed  in  this  revision.  The 
first  was  a  major  improvement  of  the  algorithm  for  classifying  which  Heppner-Maynard  high  latitude 
convection  pattern  a  given  pass  matched.  From  this  match  an  accurate  value  for  the  cross-cap  poten¬ 
tial  could  be  determined.  Thus  the  accuracy  of  the  algorithm  output  is  improved.  This  new  algorithm 
was  also  generalized  to  work  for  all  orbital  orientations,  not  just  the  dawn-dusk  orientation  of  F8. 
The  second  revision  was  a  further  refinement  of  the  endpoint  selection  procedure  to  correct  for  the 
higher  latitude  endpoints  on  the  dayside  for  the  non-dawn-dusk  oriented  orbits.  The  code  for  the 
program  was  delivered  to  Phillips  Laboratory/Gcophysics  Directorate  in  January  1993. 
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In  the  original  DMSPPOTMOD  program  each  pass  was  analyzed  to  determine  which  of  the  three 
I  leppner- Maynard  models  it  most  closely  matched.  The  algorithm  for  this  was  a  simplistic  one.  It 
assumed  that  the  satellite  travelled  only  exactly  along  the  0600—1800  hour  meridian  in  the  magnetic 
latitude/magnetic  local  time  (MLT)  coordinate  system,  then  cor.. pared  the  asymmetry  in  the  location 
of  zero  potential  point  along  that  line  in  the  three  models  to  the  asymmetry  seen  along  the  spacecraft's 
track.  The  first  failing  of  this  algorithm  comes  from  the  fact  that  the  F8  satellite  does  not  .ravel  exactly 
along  this  line  for  most  of  its  orbits  during  a  day.  Because  of  the  tilt  of  the  magnetic  dipole  axis  relative 
to  the  Earth's  rotation  axis,  the  track  of  the  F8  satellite  in  magnetic  latitudc/MLT  coordinates  “wan¬ 
ders”  between  the  nightside  and  the  06$  >- 1800  MLT  meridian  over  the  course  of  a  day  in  the  northern 
hemisphere.  In  the  southern  hemisphere,  the  ground  track  “wanders”  between  the  0600- lSOO  MI.T 
meridian  and  the  dayside  over  the  course  of  a  day.  Thus  the  difference  between  the  asymmetry  seen 
in  a  model  along  the  0600-1800  MLT  meridian  and  the  asymmetry  seen  by  the  satellite  moving  paral¬ 
lel  to  the  0600- 1800  MLT  meridian  but  15  degrees  away  from  it  towards  tne  nightside  may  be  signifi¬ 
cant.  The  second  failing  of  this  algorithm  comes  from  the  fact  that  it  docs  not  take  into  account  the 
orientation  of  the  satellite.  It  assumes  the  satellite  is  travelling  on  the  0600-  18(H)  MLT  orientation, 
but  none  of  the  DMSP  satellites  after  F8  flew  in  such  an  orientation.  Thus  the  application  of  algorithm 
to  F9,  which  flew  in  a  1000-  2200  MLT  orientation,  surely  gave  inaccurate  classifications  of  the  Hcp- 
pner-Maynard  patterns  for  a  significant  fraction  of  the  passes. 

The  offset  of  the  satellite  track  from  a  purely  0600- 1800  MLT  line  was  recognized  in  the  original 
DMSPPOTMOD  program  when  the  program  calculated  the  true  cross-cap  potential  drop  based  on 
the  observed  potential  maximum  and  minimum.  Since  it  is  unlikely  that  on  any  given  pass  the  satellite 
will  encounter  the  absolute  maximum  and  the  absolute  minimum  potential,  the  algorithm  must  some¬ 
how  scale  the  observed  maximum  and  minimum  to  get  the  cross-cap  potential.  The  original 
DMSPPOTMOD  program  attempted  this,  but  only  in  an  oversimplified  manner.  The  program  deter¬ 
mined  the  highest  magnetic  latitude  reached  by  the  satellite  on  a  given  pass,  then  used  that  to  denote 
one  of  only  four  correction  factors  (PSICORFAC)  that  were  multiplied  to  the  observed  potential  drop 
to  get  the  “true”  potential  drop.  If  the  pass  went  above  85  degrees,  then  it  assumed  that  the  observed 
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maximum  and  minimum  were  the  true  absolute  maximum  and  minimum,  so  the  correction  factor  was 
1.000.  If  the  highest  magnetic  latitude  the  satellite  reached  was  between  80  and  85  degrees,  then  the 
correction  factor  was  1.085.  If  the  highest  magnetic  latitude  the  satellite  reached  was  between  75  and 
80  degrees,  then  the  correction  factor  was  1.215.  If  the  pass  failed  to  get  above  75  degrees  magnetic 
latitude,  then  it  was  assumed  that  the  satellite  failed  to  go  through  enough  of  the  pattern  to  get  a  mean¬ 
ingful  measurement,  and  so  that  pass  was  discarded  and  the  correction  factor  was  set  to  zero.  The 
quantitative  values  of  these  correction  factors  were  based  empirically  on  the  Heppner- Maynard 
models.  The  shortcoming  of  this  algorithm  was  based  on  two  factors,  first  that  there  were  only  three 
possible  correction  factors  (along  with  the  fourth  indicating  a  discarded  pass)  instead  of  a  varying 
range  of  possible  correction  factors,  and  second,  that  it  was  based  only  on  the  dawn -dusk  orbital 
orientation  of  F8.  Any  other  orbital  orientation  required  a  completely  different  set  of  correction  fac¬ 
tors.  As  the  total  cross  -cap  potential  drop  as  well  and  the  values  of  the  absolute  maximum  and  mini¬ 
mum  potential  are  important  inputs  to  the  Magnctosphcric  Specification  and  Forecast  Model 
(MSFM)  program  being  run  by  the  Air  Weather  Service,  an  upgrade  of  this  algorithm  was  required. 
To  correct  for  this,  the  algorithm  was  redone  to  map  the  location  of  the  zero  potential  line  relative 
to  the  convection  reversal  boundary  and  the  potential  distribution  around  the  boundary  for  all  three 
models,  It  had  to  compensate  for  the  actual  track  of  the  satellite  across  the  pattern  and  use  that  to 
match  the  pass  to  one  of  the  models.  In  DMSPPOTMOD  through  DMSPDBASE3  the  model  match¬ 
ing  and  coriection  factor  analysis  were  done  in  Blocks  3F  and  3G,  respectively.  In  the  revision 
DMSPDBASF.4  the  analysis  is  combined  into  a  single  block  now  referred  to  as  Block  3F/G.  In  the 
original  program  there  were  only  four  possible  model  classifications:  1  for  model  BC,  2  for  model  A. 
3  for  model  DE.  and  4  for  unusable  pass  (probably  northward  IMF  ease).  There  are  still  only  the  three 
Heppner-Maynard  models,  but  the  new  algorithm  allows  us  to  sort  the  unusable  passes  into  several 
new  categories.  For  example,  in  the  original  algorithm,  any  pass  in  which  the  observed  potential  differ¬ 
ence  between  the  maximum  and  minimum  was  less  than  40  kV  was  classified  as  a  4.  The  rationale 
behind  this  was  that  such  a  small  drop  only  occurred  during  times  of  northward  IMF.  However  from 
the  experience  gained  by  analyzing  the  DMSP  data,  it  was  discovered  that  if  only  the  magnitude  of  the 
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potential  drop  was  known,  there  were  at  least  three  possible  cases  here:  1 )  The  satellite  passed  close 
enough  to  the  pole  so  that  it  did  likely  measure  values  close  to  the  true  potential  maximum  and  mini¬ 
mum,  and  thus  the  low  potential  drop  indicated  this  was  indeed  a  northward  IMF  case.  2)  The  satellite 
passed  close  to  the  edge  of  the  convection  reversal  boundary  which  indicated  that  the  true  potential 
drop  may  be  larger  than  40  kV,  but  the  observed  drop  was  less  than  40  k  V.  Such  passes  may  be  actually 
be  southward  IMF  cases  corresponding  to  one  of  the  three  Hcppncr- Maynard  patterns.  3)  The  satel¬ 
lite  skimmed  by  the  edge  of  the  auroral  zone  and  never  passed  through  enough  of  the  polar  region  to 
get  a  valid  measurement.  Such  passes  are  referred  to  as  “skimmers”  and  no  information  about  the 
cross-cap  potential  and  the  IMF  can  be  obtained  from  them.  Reanalyzing  the  data  using  this  revised 
algorithm  allows  us  to  differentiate  between  the  true  northward  IMF  cases  and  the  cases  of  merely 
incomplete  data.  Also,  a  few  of  the  cases  that  were  originally  classified  as  type  4  can  now  be  reclassi¬ 
fied  as  actually  one  of  the  three  Hcppncr -Maynard  patterns. 

The  algorithm  starts  by  weeding  out  the  easily  identifiable  passes.  Any  pass  with  a  potential  drop 
of  less  than  10  kV  is  identified  as  a  skimmer  and  the  model  number  (IMODNUM)  is  set  to  5.  This 
is  based  on  the  observation  that  even  during  northward  IMF  cases  where  the  satellite  passes  close  to 
the  pole,  the  observed  potential  drop  is  something  larger  than  10  kV.  Next  the  pass  is  checked  to  see 
if  cither  the  observed  maximum  or  minimum  potential  is  equal  to  zero.  If  this  is  true,  then  the  pass 
covered  a  region  of  only  positive  or  negative  potential,  and  thus  there  is  no  way  to  determine  which 
of  the  patterns  it  matches.  This  most  often  occurs  on  F9  and  FID  when  the  satellite  passes  through 
only  the  negative  potential  region.  For  this  case  the  model  number  is  set  to  6.  Next  it  examines  the 
pass  to  see  if  it  went  above  75  degree  magnetic  latitude  but  still  had  a  potential  drop  of  less  than  40 
kV.  If  this  is  true,  then  this  is  definitely  a  case  of  northward  IMF  and  the  model  number  is  set  to  4. 
If  any  of  the  above  three  cases  arc  found  to  be  true,  then  the  program  branches  to  the  end  of  the  block 
and  continue  on  with  the  rest  of  the  program.  At  this  point  all  the  passes  that  are  not  weeded  out  are 
eligible  to  be  classified  as  one  of  the  three  Hcppncr- Maynard  patterns  (Figure  5)  and  the  program 
can  now  move  on  to  the  pattern  analysis.  Some  of  the  passes  here  will  still  prove  to  be  unusable,  but 
that  will  be  taken  care  of  at  the  end  of  the  pattern  analysis  procedure. 
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(a)  (b) 


0- 

(c) 

Figure  5.  These  are  the  three  Heppner-Maynard  patterns  for  the  electrostatic  potential  dis¬ 
tribution  in  the  polar  ionosphere  (Heppner  and  Maynard ,  1987).  This  figure  shows  model  A 
( a),  model  BC  (b),  and  model  DE  (c).  The  convection  reversal  boundary  has  been  highlighted 
on  all  three  patterns.  Note  that  because  of  the  Marang  discontinuity  there  is  some  ambiguity 
as  to  the  exact  location  of  the  boundary  in  the  region  of  2100  to  0000  magnetic  local  time. 
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There  are  still  some  simplifications  in  the  revised  algorithm  as  a  full  pattern  recognition  algo¬ 
rithm  based  on  a  single  pass  is  beyond  the  scope  of  this  program  at  this  time.  The  most  significant 
simplification  is  that  we  treat  the  convection  reversal  boundary  as  a  circle.  While  the  true  boundary 
is  not  a  perfect  circle  ,  either  in  reality  or  in  any  of  the  three  Heppner- Maynard  models  (see  Figure 
5),  using  a  circle  serves  as  a  good  approximation  and  simplifies  the  analysis  procedure.  We  started 
by  defining  the  convection  reversal  boundary  and  its  center  point  for  each  of  the  three  I  leppnci  -  May¬ 
nard  models.  We  then  mapped  the  zero  crossing  line  as  a  function  of  angle  and  fraction  of  distance 
from  the  center  point  to  the  boundary  at  that  angle.  These  measurements  were  then  plotted  onto  a 
circle  which  gave  us  the  three  slightly  distorted  Heppner- Maynard  zero  potential  line  patterns  shown 
in  Figure  6.  To  determine  which  of  the  three  models  best  fits  the  observations  from  a  given  p;os,  we 
must  first  figure  out  what  the  orientation  of  the  pass  was  relative  to  the  patterns  and  where  the  pass 
crossed  the  convection  reversal  boundary.  We  assume  that  the  locations  of  the  maximum  and  mini¬ 
mum  potentials  on  a  pass  correspond  to  the  convection  reversal  boundary  and  we  convert  their  posi¬ 
tions  from  magnetic  latitude/MLT  to*  and  y  coordinates  with  the*,y  origin  set  at  the  center  of  the  mag¬ 
netic  latitude/MLT  system  (i.e. — at  90  degrees  magnetic  latitude).  We  assume  for  this  analysts  that 
the  convection  reversal  boundary  is  a  perfect  circle  that  is  centered  near,  but  not  necessarily  at,  the 
coordinate  origin.  Since  there  are  only  two  points  known  here  (the  maximum  and  minimum  locations) 
there  are  an  infinite  number  of  circles  which  contain  both  these  points  on  their  circumference.  To 
choose  a  single  circle  we  make  a  second  simplifying  assumption:  that  the  radius  of  the  circle  is  equal 
to  the  average  of  the  distances  between  each  extrema  point  and  the  coordinate  origin.  Again,  this  is 
not  exact,  but  in  most  eases  it  gives  a  good  approximation  of  the  true  location  of  the  convection  reversal 
boundary.  With  the  radius  of  the  circle  given  and  the  locations  of  the  maximum  and  minimum  known, 
we  have  narrowed  down  the  location  of  the  center  of  the  circle  to  only  two  points.  The  algorithm  calcu¬ 
lates  the  locations  of  those  two  points  in  x,y  coordinates,  then  selects  the  one  closest  to  the  coordinate 
origin  to  be  the  center  of  the  circle. 

With  the  radius  and  location  of  the  center  of  the  model  convection  reversal  boundary  determined, 
the  program  converts  the  locations  of  the  observed  maximum,  minimum,  and  zero  potential  into  the 
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Figure  6.  These  patterns  show  the  simplified  convection  reversal  boundary  and  the  zero  poten¬ 
tial  line  used  in  the  program  to  identify  the  Heppner- Maynard  model  A  (a),  model  BC  (b)t 
and  model  DE  (c).  These  patterns  have  been  distorted  slightly  relative  to  the  true  patterns 
shown  in  figure  5  by  forcing  the  boundary  to  be  a  true  circle.  The  asterisks  on  the  boundary 
denote  the  potential  distribution  along  the  boundary  in  steps  of  4  kV. 
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new  coordinate  system  where  the  center  of  the  circle  is  the  origin.  It  then  calls  three  subroutines  ( MO- 
DELA,  MODELBC,  and  MODELDE)  to  calculate  the  location  of  the  zero  potential  point  along  the 
satellite's  track  and  the  correction  factors  for  the  maximum  and  minimum  potentials  for  each  of  the 
three  Heppner-Maynard  models.  Each  of  these  subroutines  has  only  the  data  arrays  for  that  specific 
model  and  that  it  then  calls  a  further  subroutine  named  MODEL.  to  do  the  actual  work.  This  allows 
for  future  models  to  be  easily  inserted  into  the  program  by  simply  adding  another  subroutine  (say  MO- 
DELFG)  that  contains  only  the  data  arrays  to  specify  the  potential  distribution  around  the  convection 
reversal  boundary  and  the  location  of  the  zero  potential  line  for  that  model.  The  location  of  the  ob¬ 
served  zero  potential  point  along  the  track  is  used  to  calculate  the  distance  between  it  and  the  location 
of  the  maximum.  That  distance  is  divided  by  the  calculated  length  of  the  chord  between  the  locations 
of  the  maximum  and  minimum  potential  to  get  the  fractional  distance  called  FRACORIG.  The  pro¬ 
gram  then  takes  the  returned  locations  of  the  zero  potential  point  for  each  of  the  three  Heppner- 
Maynard  models,  then  calculates  their  fractional  distances  along  the  chord  of  the  satellite  track  (FRA- 
CA,  FRACBC,  and  FRACDE).  Finally  the  program  chooses  the  model  with  the  fractional  distance 
that  is  closest  to  the  observed  fractional  distance  as  the  model  that  best  fits  this  pass.  The  program 
sets  the  model  number  to  1  for  pattern  BC,  2  for  pattern  A,  and  3  for  pattern  DE,  which  arc  the  same 
numbers  used  in  the  original  DMSPPOTMOD  program.  At  this  point  the  program  can  weed  out 
some  of  the  ambiguous  or  distorted  passes.  Examination  of  the  patterns  in  Figure  6  show  that  for  al¬ 
most  every  possible  orientation  of  the  satellite  track,  the  smallest  fractional  distance  from  the  ob¬ 
served  maximum  to  the  observed  zero  potential  point  always  belongs  to  model  BC  and  the  largest  al¬ 
ways  belongs  to  model  DE.  Only  in  two  regions  docs  this  rule  fail.  On  the  dayside  between  09(X)  and 
1300  MLT  and  below  60  degrees  magnetic  latitude  the  location  of  the  zero  lines  in  model  A  and  model 
DE  trade  places.  On  the  nightsidc  near  0000  MLT  and  at  magnetic  latitudes  less  than  60  degrees,  zero 
potential  lines  get  so  close  together  that  it  become  impossible  to  determine  unambiguously  which  of 
the  three  models  is  the  best  fit.  If  the  track  of  the  satellite  is  such  that  it  intersects  cither  of  these  re¬ 
gions  then  the  returned  values  for  the  fractional  distances  from  the  three  models  fail  to  meet  the 
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FRACBC  <  FRACA  <  FRACDF  condition.  For  those  cases  no  determination  can  he  made  unambig¬ 
uously,  so  those  cases  are  ruled  as  unusable  and  the  model  number  is  set  to  4. 


When  the  MODELA.  MODELBC,  and  MODELDE  subroutines  arc  called  they  each  return  two 
correction  values.  These  corrections  are  the  numbers  which  arc  multiplied  to  the  observed  magni¬ 
tudes  of  the  potential  maximum  and  minimum  to  get  the  absolute  potential  maximum  and  minimum 
for  that  model.  These  correction  factors  are  calculated  by  taking  the  potential  distribution  around 
the  convection  boundary  for  the  given  model  and  comparing  the  potential  at  the  location  of  the  ob¬ 
served  maximum  (minimum)  to  the  absolute  maximum  (minimum)  for  that  model.  When  the  pro¬ 
gram  chooses  the  best  fit  model,  then  the  correction  factors  for  that  model  are  set  to  the  variables 
CORRECTMAX  and  CORRECTMIN.  For  the  northward  IMF,  skimmer,  and  other  unusable 
passes,  the  values  of  CORRECTMAX  and  CORRECTMIN  are  set  to  zero.  The  greater  the  distance 
betw  een  the  location  of  the  observed  maximum  (minimum)  and  the  location  of  the  absolute  maximum 
(minimum)  in  the  model,  the  larger  the  correction  is  going  to  be  and  the  less  certainty  such  a  correction 
carries.  For  example,  an  observed  maximum  of  65  kV  near  0700  MLT  with  a  correction  factor  of  1.14 
is  obviously  more  certain  than  an  observed  maximum  of  7  kV  at  1030  MLT  with  a  correction  value  of 
10.57.  While  we  cannot  at  this  time  give  any  quantitative  reliability  to  these  correction  factors  (other 
than  the  general  caveat  that  the  larger  the  correction,  the  more  wary  the  investigator  should  be  about 
using  it)  we  have  set  a  value  of  10  for  the  cutoff  for  the  correction  factors.  Any  pass  that  ends  up  with 
a  correction  factor  of  10  or  more  is  ruled  unusable  and  the  model  number  is  set  to  7  and  the  values 
of  CORRECTMAX  and  CORRECTMIN  arc  set  to  zero.  Also,  there  arc  a  few  cases  where  the  real 
pattern  is  so  distorted  that  the  observed  maximum  (minimum)  occurs  in  a  region  of  the  boundary 
which  the  model  designates  as  negative  (positive)  and  thus  returns  a  negative  value  for  the  correction 
factor.  Such  a  case  is  too  distorted  to  be  used  and  is  also  classified  with  a  model  number  of  7.  This 
gives  us  a  total  of  eight  possible  categories  for  classifying  each  pass.  The  eight  categories  are  summa¬ 
rized  in  the  table  below. 

Table  of  model  numbers  for  DMSPDBASE4 
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O*northward  IMF  (delta  PSI  <  40  kV,  highest  magnetic  latitude  point  >  75  degrees) 

1  HM  model  BC  \ 

2  HM  model  A  I — determined  from  model  fitting  subroutines 

3  HM  model  DE  / 

4  unusable  (zero  occurred  in  far  nightside  or  else  pattern  too  distorted  to  classify  unambiguously) 

5  skimmer  (unusable — delta  PSI  <  10  kV) 

6  unusable  (either  observed  maximum  or  minimum  was  zero) 

7  unusable  (one  of  the  corrections  greater  than  10  or  negative,  pattern  likely  too  distorted  to  classify 

unambiguously) 


*Note  that  for  null  passes  the  value  for  the  model  is  also  set  to  zero.  However  these  cases  are 
easily  distinguished  from  the  northward  IMF  cases  in  that  null  passes  have  zero  values  for  all  the 
parameters. 

We  have  compared  the  model  classification  results  of  two  sets  of  ten -day  periods  of  F8  and  Fit) 
data  using  the  original  DMSPPOTMOD  through  DMSPDBASE3  algorithm  to  the  new  model  classifi¬ 
cation  results  using  DMSPDBASE4  and  summarized  the  results  in  the  tables  below.  Note  that  for 
F8,  there  is  not  much  change  between  the  number  of  passes  classified  as  a  Hcppncr-  Maynard  pattern 
BC  (model  1)  and  pattern  DE  (model  3)  between  the  two  algorithms.  About  half  of  the  passes  origi¬ 
nally  classified  as  pattern  A  (model  2)  were  reclassified  as  being  pattern  BC.  As  the  original  method 
assumed  an  F8pass  on  a  nearly  dawn -dusk  track,  the  fact  that  the  old  algorithm  performed  fairly  well 
on  F8  data  is  not  surprising. 


rx  i- 

10  January  1989 

orig 

#  of  passes 

new  0 

new  1 

new  2 

new  3 

new  4 

new  5 

new  6 

new  7 

null 

1 

68 

3 

51 

2 

1 

5 

0 

2 

4 

— 

2 

33 

3 

14 

12 

1 

2 

0 

0 

1 

— 

3 

4 

0 

0 

1 

3 

0 

0 

0 

0 

— 

4 

157 

69 

10 

1 

0 

11 

43 

21 

4 

— 

null 

13 

3 

1 

0 

0 

0 

1 

0 

0 

8 

total 

277 

78 

76 

16 

5 

IX 

44 

23 

9 

8 

F8  1- 

■10  April  1989 

orig 

#  of  passes 

new  0 

new  1 

new  2 

new  3 

new  4 

new  5 

new'  6 

new  7 

null 

1 

2 


104 

76 


81 

32 


4 

36 


0 

1 


10 

3 


0 

0 


6 

0 


2 

2 


2 


29 


total 


43 


5 


0 

28 

0 


0 


6 

87 

3 


0 

32 

0 


0 

7 

2 


-> 

4- 

0 


0 

9 

0 


0 

It) 

0 


t) 

2 

0 


276 


35  118 


22 


28 


16 


The  next  two  tables  show  the  results  of  the  new  algorithm  on  two  sets  of  F10  data.  Since  Fit)  is 
in  a  08(X)-2000  MLT  orbit  the  original  algorithm  did  not  perform  as  well  as  it  did  for  F8.  After  using 
the  new  analysis  only  about  45f «  of  the  passes  originally  classified  as  pattern  BC  (model  1)  remain 
as  pattern  BC,  while  about  25%  arc  reclassified  as  pattern  A  (model  2)  and  the  remainder  fall  into 
one  of  the  non  -Heppner-  Maynard  categories.  There  was  less  change  seen  in  the  pattern  A  (model 
2)  passes.  Just  over  60%  of  the  passes  originally  classified  as  pattern  A  remained  in  that  category  after 
the  new  analysis,  while  about  8%  were  reclassified  as  pattern  BC,  and  about  10%  were  reclassified 
as  pattern  DE  (model  3).  Overall  only  around  one-third  of  the  F10  passes  fall  into  one  of  the  three 
Heppner-Maynard  patterns  compared  to  nearly  one-half  of  the  F8  passes,  while  the  percentage  of 
“unsuitable  passes”  (model  numbers  4-7)  increased  from  30%  for  F8  to  38%  for  F10.  This  is  likely 
a  result  of  the  different  orbital  orientation  of  the  F10  spacecraft  relative  to  F8. 


F10  1 

-If)  January  1991 

orig 

#  of  passes 

new  0 

new  1 

new  2 

new  3 

new  4 

new  5 

new  6 

new  7 

null 

1 

33 

3 

15 

9 

0 

3 

0 

1 

2 

— 

2 

24 

0 

1 

15 

5 

0 

0 

0 

3 

— 

3 

13 

0 

0 

2 

10 

1 

0 

0 

0 

4 

153 

77 

3 

6 

1 

19 

38 

7 

2 

— 

null 

26 

4 

3 

2 

3 

3 

4 

1 

2 

4 

total 

249 

84 

22 

34 

19 

26 

42 

9 

9 

4 

F10  1-10  April  1991 

orig  #  of  passes 

new  0  new  1 

new  2 

new  3 

new  4 

new  5  new  6 

new  7  null 

1  76 

2  34 

18 

0 

3 

0  17 

2 
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2  24  0  3  14  0  1  0  0  (. 

3  9  0  0  15  2  0  0  1 

4  93  34  I  !  1  2X  3  24  1 

null  42  7  5  5  1  6  3  7  0  X 

Total  244  43  43  39  7  40  6  4X  5Ti  x 


An  additional  benefit  of  the  new  algorithm  can  be  seen  in  the  fact  that  a  large  number  of  the  passes 
originally  classified  as  “null  passes”  under  the  old  algorithm  turn  out  to  be  analyzablc  passes  under 
the  new  algorithm,  although  about  half  of  them  still  end  up  classified  in  one  of  the  “unsuitable  passes" 
categories.  This  is  particularly  dramatic  in  the  two  F10  tables  where  over  80%  of  the  "null  passes" 
under  the  old  algorithm  become  good  passes  under  the  new  algorithm.  This  raises  the  percentage  of 
analyzable  passes  for  the  two  ten-day  periods  of  FID  from  the  80%  range  to  the  90%  range 

The  requirement  that  the  potential  correction  factors  be  saved  to  the  database  along  with  the  re¬ 
doing  of  the  model  number  classification  scheme  meant  that  we  had  to  reformat  the  shortfilc  To  keep 
the  size  of  the  128  character  record  per  pass  constant  we  had  to  delete  something  to  make  room  for 
the  correction  factors.  The  value  of  the  average  of  the  potential  maximum  and  minimum  (POT- 
MEANSF)  and  the  location  where  it  occurred  (POTMEANMLT  and  IPOTMEANMLAT)  had  not 
proven  to  be  a  useful  parameter.  Furthermore,  since  the  entire  set  of  potentials  for  a  given  pass  wt  c 
being  saved  in  the  longfilc,  these  values  could  be  easily  recovered  if  the  need  ever  arose.  These  thr.e 
values  were  replaced  with  the  two  correction  values  (CORRECTMAX  and  CORRECTMIN)  and  the 
highest  magnetic  latitude  the  satellite  reached  on  this  pass  (MLATHIGH).  As  with  the  other  magnetic 
latitudes  in  the  shortfilc,  the  value  saved  is  the  latitude  times  ten  in  order  to  remove  the  need  for  the 
decimal  point.  This  variable  will  enable  investigators  to  easily  ascertain  the  degree  to  which  any  given 
pass  penetrated  the  polar  cap  region. 

IQUALFLAG  is  now  a  only  2-digit  integer,  but  still  formatted  as  14  as  in  the  original  shortfilc. 
The  hundreds  place  was  the  IPSICORFAC  (PSf  CORrection  FACtor)  variable,  which  no  longer  exists 
since  it  is  now  superceded  by  the  new  CORRECTMAX  and  CORRECTMIN  variables  .  The  tens 
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place  is  still  IQFLAG  (Quality  FLAG)  which  is  unchanged  from  the  original  shortfile.  The  final  digit 
is  the  model  number  (IMODNUM)  classification  discussed  above  and  now  ranges  from  0  to  7. 


The  two  data  file  formats  for  Urn  shnrtfile  are  given  below  (the  *  indicates  a  change  between  the 


two): 


OLD  SHORTFILE  FORMAT  NEW  SHORTFILE  FORMAT 


item 

format 

item 

format 

SFINDEX 

A12 

SFINDF.X 

A12 

PSIMAXSF 

F5.1 

PSIMAXSF 

F5.1 

SCMLTMAX 

F5.1 

SCMLTMAX 

F5.1 

1NVLATMAX 

15 

INVLATMAX 

15 

I  ML  ATM  AX 

15 

IMLATMAX 

15 

PSIMINSF 

F7.1 

PSIMINSF 

F7.1 

SCMLTMIN 

F5.1 

SCMLTMIN 

F5.1 

INVLATMIN 

15 

INVLATMIN 

J5 

1MLATMIN 

15 

IMLATMIN 

15 

IQl'ALFLAG 

14 

* 

IQUALFLAG 

14 

ZEROMLT 

F5.1 

ZERGMLT 

F5.I 

IZEROMLAT 

15 

IZEROMLAT 

15 

PLMIDPOTSF 

F7.I 

PLMIDPOTSF 

F7.1 

PLMIDMLT 

F5.I 

PLMIDMLT 

F5.1 

1PLMIDMLAT 

15 

IPLMIDMLAT 

15 

POTMEANSF 

F7.1 

* 

CORRECTMAX  F6.3 

POTMEANMLT 

F5. 1 

* 

CORRECTMIN 

F6.3 

I POTM  E  A  N  ML  AT 

15 

* 

MLATHIGH 

15 

KPSHORT 

13 

KPSHORT 

13 

IAE1NDEX 

14 

IAE1NDEX 

14 

BXSHORT 

F5.1 

BXSHORT 

F5.1 

BYSHORT 

F5.1 

BYSHORT 

F5.1 

BZSHORT 

F5. 1 

BZSHORT 

F5.1 

1POTOFF 

14 

IPOTOFF 

14 

Total  24  columns  128  characters  plus  oneOxOa  code  (a  linefeed  at  the  beginning  of  each  line),  form¬ 
ing  129  bytes/record. 


The  second  major  revision  in  DMSPDBASE4  is  a  further  refinement  of  the  endpoint  selection 
algorithm.  As  was  explained  above,  the  revised  algorithm  in  DMSPDBASF.2  began  looking  for  each 
endpoint  by  examining  the  data  at  65,  60,  55,  and  50  degrees  magnetic  latitude  and  if  more  than  one 
was  acceptable,  it  choose  the  one  that  was  the  lowest  in  latitude.  This  choice  was  based  on  the  assump- 
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tion  that  if  a  higher  latitude  point  was  also  acceptable,  then  the  potential  drop  between  the  higher  and 
lower  latitude  points  was  zero,  and  that  extending  the  “tail”  of  the  potential  curve  at  either  end  would 
only  result  in  a  flat  line  at  essentially  zero  potential.  For  most  eases  this  assumption  held,  as  can  be 
seen  in  Figure  7  where  the  potential  is  essentially  zero  from  50  degrees  to  60  degrees  magnetic  latitude 
on  the  dawnsidc  (left  side  of  graph)  and  from  around  57  degrees  to  50  degrees  on  the  duskside.  The 
preference  of  50  degrees  magnetic  latitude  as  a  starting  and  ending  latitude  was  historical  in  its  basis. 
Fifty  degrees  as  a  nominal  and  safe  choice  for  the  zero  potential  was  used  throughout  all  the  earlier 
work  on  NASA’s  Dynamics  Explorer.  Most  of  the  concern  had  been  on  making  the  algorithm  robust 
enough  to  handle  eases  when  the  polar  cap  boundary  expanded  during  storm  periods  where  the  start¬ 
ing  and  stopping  endpoints  had  to  be  reset  cquatorward  of  the  50  degree  line.  However  as  the  analysis 
progressed  there  were  some  occasions  where  a  combination  of  a  small  polar  cap  boundary  combined 
with  differences  in  the  flow  baseline  on  either  side  of  the  pass  would  result  in  a  false  potential  trough 
or  hump  appearing  cquatorward  of  the  auroral  region.  Figure  8  shows  an  example  of  this  on  the  dusk- 
side  (right  side)  of  the  potential  curve.  The  slight  difference  in  the  baselines  in  the  horizontal  flow 
on  either  side  of  the  polar  region  would  not  make  very  much  difference  if  the  endpoints  were  just 
cquatorward  of  the  auroral  regions.  But  the  slight  offset  of  the  flow  during  the  six  minute  period  from 
minute  46  through  minute  52  (roughly  65  degrees  magnetic  latitude  to  50  degrees)  results  in  a  shallow- 
basin  of  negative  potential  in  the  corrected  potential  curve.  Since  this  occurs  on  the  dawnsidc  of  a  track 
that  is  almost  exactly  along  the  06(H)- 1800  line,  this  portion  of  the  potential  should  normally  be  either 
positive  or  zero,  and  thus  this  negative  potential  basin  is  clearly  an  artifact  of  the  processing  procedure. 
While  this  occurred  infrequently  in  the  F8  data,  it  appeared  more  otten  in  the  F9  and  t- 10  data.  The 
region  of  zero  potential  reaches  to  higher  latitudes  on  the  dayside  than  on  the  nightside  or  on  the 
dawn -dusk  line,  thus  the  tilt  of  these  satellites’  orbits  away  from  the  dawn -dusk  orientation  of  F8 
causes  them  to  observe  zero  flows  (i.c.  zero  potential)  on  the  dayside  at  higher  latitudes  than  ob¬ 
served  by  F8.  These  artifacts  interfere  with  the  accuracy  of  electric  field  modelling  routine  in  the 
MSFM,  and  this  was  brought  to  our  attention  by  the  personnel  from  Rice  University  at  the  Quarterly- 
Review  Meeting  held  in  Colorado  Springs  in  November  1992  (see  Freeman  et  al„  1992).  The  MSFM 
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Figure  7.  I  his  is  u  plot  ot  a  pass  where  thcac  is  a  laiee  amount  of  near-  zero  flow  between  50  and 
60  decrees  on  both  sides  of  the  polar  region.  Note  that  these  tails"  nave  little  effect  on  the  potential 
curve  calculated  and  shown  in  the  bottom  panel 
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Figure  8.  This  pass  is  similar  to  the  one  in  figure  7.  except  that  the  flow  equatorward  of  the  auroral 
region  is  does  not  go  to  zero.  Since  there  is  a  difference  in  the  baseline  determined  by  the  flow  on 
the  duskside  (left)  and  the  dawnside  (right),  the  program  averages  the  difference  and  comes  up  with 
a  baseline  that  gives  non-zero  flows  on  either  side.  This  small,  but  non-zero  flow  in  these  “tails” 
results  in  an  artificial  negative  potential  region  on  the  dawnside  between  the  auroral  region  and  50 
degrees. 


attempts  to  fit  a  potential  curve  to  the  erroneous  portion  of  a  potential  curve  such  as  shown  in  Figure 
8  which  obviously  degrades  the  ability  of  the  MSFM  to  produce  quality  results.  At  their  request  we 
began  to  redo  the  algorithm  to  take  care  of  these  artifacts.  Initially  there  was  some  concern  that  any 
change  that  would  push  the  endpoint  towards  higher  latitudes  for  these  problem  passes,  might  inad¬ 
vertently  degrade  the  quality  of  the  analysis  of  the  nominal  passes.  By  pushing  the  endpoints  of  nomi¬ 
nal  passes  to  higher  latitudes  this  might  change  the  baselines  which  would  change  the  calculated  po¬ 
tentials,  or  it  might  even  accidently  push  the  endpoints  into  the  polar  region  and  miss  the  auroral 
portion  of  the  pass  altogether. 

Fortunately,  the  fix  turned  out  to  be  quite  simple  and  caused  little  or  no  perturbation  of  the  data 
for  the  nominal  passes.  The  revised  algorithm  gets  rid  of  the  50  degree  preference  by  simply  starting 
the  check  procedure  of  the  data  at  65  degrees  magnetic  latitude  and  proceeding  cquatorward  in  five 
degree  ctcps  until  either  an  acceptable  endpoint  or  20  degrees  is  reached.  This  guarantees  that  the 
highest  latitude  point  that  is  suitable  is  used  as  the  endpoint  for  that  end  of  the  pass.  Figure  9  shows 
the  new  potential  for  the  pass  shown  in  Figure  8  using  this  new  algorithm  to  choose  the  endpoints. 
This  time  the  pass  starts  at  60  degrees  magnetic  latitude  and  stops  at  65  degrees.  (Note  that  in  this 
report  and  in  the  program  we  use  only  the  absolute  value  of  magnetic  latitude  for  our  analysis  and 
descriptions.  The  negative  sign  on  the  magnetic  latitude  shown  in  the  figure  above  merely  indicates 
that  the  pass  occurred  in  the  southern  hemisphere.)  Note  that  this  gets  rid  of  the  false  negative  poten¬ 
tial  region  on  the  dawnsidc  and  resets  the  maximum  and  minimum  potential  values,  probably  closer 
to  the  true  values  than  was  the  ease  in  the  first  analysis.  Figure  10  presents  the  same  pass  as  was  shown 
in  Figure  7  after  using  the  new  algorithm.  This  results  in  only  clipping  the  zero  potential  regions  on 
either  end.  The  corrected  maximum  and  minimum  potentials  arc  changed  by  less  than  2%  from  their 
values  in  the  old  analysis.  While  this  new  version  decreases  the  likelihood  of  faise  potentials  appearing 
in  the  data  because  of  a  bad  choice  of  endpoints,  it  should  he  strongly  emphasized  that  this  routine 
is  still  not  absolutely  100%  foolproof.  The  level  of  confidence  in  all  of  this  analysis  must  continually 
be  reevaluated  in  light  of  the  changing  conditions  of  the  data  that  the  program  faces.  As  we  move 
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Figure  9.  This  is  the  same  pass  as  in  figure  X  except  that  the  endpoints  have  been  moved  poleward 


to  avoid  including  the  non-zero  “tails”  in  the  flow  in  the  potential  calculations.  Note  that  the  artifi¬ 


cial  negative  potential  region  has  disappeared.  Also  the  magnitudes  of  the  maximum  and  minimum 


and  location  of  the  zero  point  have  changed  from  figure  X. 
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Figure  10.  This  is  the  same  pass  as  in  figure  7  except  that  the  endpoints  have  been  moved  poleward 


to  avoid  including  the  near  zero  "tails”  in  the  potential  calculation.  Note  that  removing  the  "tails' 
has  not  significantly  changed  the  locations  or  magnitudes  of  the  maximum  or  minimum. 


towards  solar  minimum,  the  ionospheric  conditions  will  likely  change  and  this  may  be  reflected  in  the 
quality  of  the  analysis,  which  in  turn  may  require  some  further  upgrades  be  made  to  DMSPDBASE4. 


4.  OTHER  ANALYSIS  EFFORTS 

While  most  of  the  effort  performed  on  this  contract  has  been  on  upgrading  and  maintaining  the 
analysis  programs  we  have  developed,  we  have  devoted  a  fair  amount  of  time  to  efforts  using  the  data 
generated  for  the  DMSP  database.  Our  major  focus  has  been  on  characterizing  the  convection  pat¬ 
terns  observed  during  times  of  steady  IMF  conditions  when  Bz  is  southward.  With  the  continuous  data 
stream  from  DMSP,  we  can  use  take  extended  times  of  steady  IMF  conditions  to  compare  the  patterns 
seen  simultaneously  in  both  polar  region  using  multiple  satellites.  Generally  we  have  found  that  the 
patterns  seen  are  consistent  with  the  Hcppner- Maynard  patterns,  but  we  have  been  surprised  to  find 
that  there  is  roughly  a  10%  difference  between  the  total  cross -cap  potential  drops  observed  in  the 
two  hemispheres.  We  are  in  the  final  stages  of  preparing  a  paper  on  our  results  of  this  study  for  submis¬ 
sion  to  the  Journal  of  Geophysical  Research.  In  conjunction  with  this  work  we  have  developed  a  three- 
dimensional  presentation  of  the  data  which  allows  the  user  to  examine  the  magnitude  and  shape  of 
the  potential  curve  at  the  same  time  as  the  location  of  the  satellite’s  track  in  the  magnetic  latitude/mag- 
nctic  local  time  coordinate  system  is  shown  (Figure  1 1).  This  enables  us  to  clearly  present  a  larger 
set  of  data  in  a  single  figure  without  overwhelming  the  audience  or  reader.  A  short  animation  of  the 
data  using  this  technique  was  submitted  on  videotape  to  Phillips  Lab  in  spring  1991. 

Judy  Cumnock,  one  of  the  graduate  students  here  at  the  University  of  Texas  at  Dallas,  has  been 
concentrating  on  identifying  flow  patterns  observed  during  times  of  extended  northward  IMF  condi¬ 
tions.  Her  work  studying  the  pattern  during  the  January  14,  1988  period  showed  evidence  of  a  two- 
ccll  pattern  metamorphosizing  into  a  four-cell  pattern  as  the  IMF  orientation  rotates  from  predomi¬ 
nantly  By  to  predominantly  +BZ.  This  work  has  been  presented  at  an  American  Geophysical  Union 
meeting  and  has  been  published  in  the  Journal  of  Geophysical  Research  (Cumnock  et  al„  1992).  She 
is  expanding  this  work  into  a  doctoral  thesis  which  will  examine  more  cases  of  convection  patterns  dur¬ 
ing  northward  IMF  conditions  for  the  1988-1990  period.  We  currently  have  two  other  graduate  stu¬ 
dents  working  on  analysis  of  the  DMSP  data  for  their  doctoral  work.  Chris  Keating  has  worked  on 
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Figure  1 1 .  This  is  a  three  dimensional  representation  of  a  polar  pass.  The  viewpoint  is  from  the 
nightside  looking  towards  the  dayside.  The  ground  track  is  represented  by  the  dashed  line  on  the 
dial.  The  potential  is  represented  by  the  curve  above  (positivfrpotential)  and  below  (negative  poten¬ 
tial)  the  dashed  line.  The  locations  and  magnitudes  of  the  potential  maximum  and  minimum  are 
shown  by  the  two  vertical  lines  connecting  the  potential  curve  and  the  ground  track.  In  front  of  the 
noon-midnight  line  are  two  glyphs.  The  one  on  the  left  is  a  sliding  bar  that  shows*he~Kjrvalue 
during  the  pass.  The  other  glyph  shows  the  orientation  and  magnitude  of  the  IMF  in  the  y-z  plane. 
Note  that  this  is  oriented  such  as  it  would  appear  to  an  observer  looking  towards  the  sun  from  the 
Earth.  Thus  y  is  positive  to  the  left  and  negative  to  the  right.  Here  we  show  a  case  where  Bz  and 
By  are  both  negative. 
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identifying  periods  when  'he  orbits  of  DMSP-F8  and  F9  were  in  phase  with  each  other,  then  using 
the  data  from  those  passes  to  characterize  the  shape  of  the  convection  reversal  boundary.  Keith  West 
just  started  this  past  year  and  is  looking  at  equatorial  data  from  DMSP,  a  portion  of  the  database  that 
we  have  generally  ignored  up  to  this  point.  Both  Cumnock’s  and  West  s  work  on  DMSP  analysis  have 
been  funded  by  the  Texas  Space  Grant  Consortium. 

A  major  portion  of  our  effort  in  data  analysis  went  into  producing  and  providing  usable  DMSP 
data  for  several  other  investigators  and  other  institutions.  Feedback  from  these  other  groups  has  been 
invaluable  in  honing  the  computer  analysis  tools.  As  the  software  has  matured,  fulfilling  these  re¬ 
quests  has  become  more  routine  and  automated.  However  some  of  them  still  require  special  handling 
to  produce,  such  as  for  cases  where  the  user  wants  only  data  that  was  acquired  only  over  certain  geo¬ 
graphic  regions.  Wc  have  provided  data  to  the  following  individuals  or  groups: 

*  Major  Delores  Knipp  at  the  US  Air  Force  Academy.  She  has  incorporated  the  measured  potential 
from  DMSP  into  the  Assimilative  Mapping  of  Ionospheric  Electrodynamics  (AMIE)  routine  to  help 
in  modelling  the  polar  electric  field. 

*  Drs.  Richard  Wolf  and  Robert  Spiro  at  Rice  University.  They  have  developed  the  MSFM  package 
for  use  by  the  Air  Weather  Service  and  much  of  our  work  described  above  on  DMSPDBASE  was  moti¬ 
vated  by  the  needs  of  the  MSFM  program. 

*  Dr.  John  Freeman  and  Akira  Nagai  at  Rice  University.  They  have  been  working  on  developing  a 
neural  network  which  will  predict  the  future  convection  pattern  in  the  polar  ionosphere  based  on  the 
previous  several  hours’  observations.  This  work  uses  the  shortfilc  output  from  DMSPDBASE  as  a 
testbed. 

*Dr.  Mark  Loranc  at  Marshall  Space  Flight  Center  (formerly  at  the  University  of  Western  Ontario). 
He  has  been  comparing  the  flow  velocities  observed  by  DMSP  when  it  passes  over  the  BARS  radar 
site  in  Canada  with  the  measured  ion  flows  seen  by  the  radar. 

*Dr.  Ennio  Sanchez  at  Johns  Hopkins  Applied  Physics  Laboratory.  He  has  used  DMSP-  F9  ion  flow- 
data  to  compare  to  and  calibrate  data  from  the  PACE  radars  during  substorms. 

*  Dr.  David  Knudsen  at  Max— Planck -Institut  fur  extraterr  Physik.  He  compares  the  DMSP  data 
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with  the  radar  data  from  the  EISCAT  radar. 

*Dr.  Geoff  Crowley  at  Lowell  Univeristy.  He  compares  the  DMSP  flow  with  radar  data  in  the  cusp 
region. 

*Drs.  Barbara  Emery  and  Gang  Lu  at  the  National  Center  for  Atmospheric  Studies  (NCAR).  They 
have  used  the  ion  flow  from  several  of  the  DMSPs  to  do  ionospheric  modelling. 

*  Dr.  Mervyn  Freeman  at  the  British  Antarctic  Survey.  He  and  his  colleagues  used  DMSP  flow  data 
to  model  the  ionospheric  polar  convection  pattern  during  specific  times  of  northward  IMF. 

Much  of  this  work,  both  here  at  UT  Dallas  and  in  collaboration  with  other  groups,  has  been  pres¬ 
ented  at  meetings  of  the  American  Geophysical  Union.  Below  is  a  list  of  the  AGU  presentations  over 
the  past  three  years  which  have  dealt  with  DMSP  data  analysis  covered  in  this  contract.  An  asterisk 
indicates  a  presentation  given  by  a  member  of  the  research  group  here.  The  final  letter  and  number 
code  after  the  title  gives  the  AGU  code  number  for  that  presentation. 

AGU  Spring  1990 

none 

AGU  Fall  1990 

‘Hairston,  M.  R.,  R.  A.  Heclis,  and  F  J.  Rich,  Characterization  of  Ionospheric  Flow  Patterns  in 
Northern  and  Southern  Polar  Regions  Using  DMSP  Data  SM12B-2 
Knipp.  D.  J.,  B.  A.  Emery,  A.  D.  Richmond,  M.  R.  Hairston,  R.  A.  Hcclis,  and  F.  J.  Rich,  Mapping 
Ionospheric  Convection  with  Satellite  Ion  Drift  Measurements  SM12B-  1 

AGU  Spring  1991 

‘Hairston,  M.  R.,  R.  A.  Heclis,  and  F.  J.  Rich,  Comparison  of  Simultaneous  Ionospheric  Flow  Pat¬ 
terns  in  Norncrn  and  Southern  Polar  Regions  Using  DMSP  Data  SA51B-4 
Rich,  F.  J.,  W.  J.  Dcnig,  M.  R.  Hairston,  and  R.  A.  Hcclis,  A  Search  for  the  Signature  of  Flux  Transfer 
Events  Near  the  Daysidc  Ionospheric  Cusp  SM42A-7 


“Hairston,  M.  R.,  R.  A.  Hcclis,  and  F.  J.  Rich,  Ionospheric  Polar  Cap  Potential  Distribution  for  March 
1 989  -  -  The  Movie  S  A 1 2  A  -  7 

“Cumnock,  J.  A.,  R.  A.  Hcclis,  and  M.  R.  Hairston,  Ionospheric  Ob>crvations  of  Sunward  Convection 
at  Very  High  Latitudes  SA12A-9 

AGE  Spring  1992 

Nagai,  A.,  J.  W.  Freeman  Jr.,  P.  H.  Reiff,  R.  A.  Wolf,  F.  J.  Rich,  and  \1.  R.  Hairston,  Short-term  Fore¬ 
cast  of  Magnctosphcric  Parameters  by  a  Neural  Network  System  SM32A-3 

AGU  Fall  1992 

F.mery,  B.  A.,  A.D.  Richmond,  D.  J.  Knipp,  M.  R.  Hairston,  and  W.  F.  Denig,  Polar  Ionospheric  Elec¬ 
trodynamics  and  Energy  Inputs  in  the  Northern  and  Southern  Hemispheres  for  a  Strong  I  nterpla- 
nctaty  Magnetic  Field  SA21 B-2 

Rich,  F.  J,  M.  S.  Gussenhoven,  and  M.  R.  Hairston,  DMSP  Observation  of  Precipitating  Particles  and 
Electric  Field  During  the  GEM  Campaign  SM41B-4 

Boyle,  C.  B,  P.  H.  Reiff,  A.  Nagai,  M.  R.  Hairston,  and  R.  A.  Hcclis,  DMSP  Driven  Estimates  of  Polar 
Cap  Potential  SM51A-1 

Much  of  this  work  has  been  presented  as  research  papers  which  arc  in  various  states  of  acceptance 

in  the  refereed  literature.  The  list  of  papers  is  presented  below: 

Papers  published 

Cumnock,  J.  A.,  R.  A.  Hcclis,  and  M.  R.  Hairston,  Response  of  the  ionospheric  convection  pattern 
to  a  rotation  of  the  interplanetary  magnetic  field  on  January  14,  1988,/  Geophys.  Res.,  97, 
19449-19460,  1992. 

Papers  accepted  and  awaiting  publication 

Freeman,  M.  P.,  C.  J.  Farrugia,  L.  F.  Burlaga,  R.  P.  Lepping,  M.  E.  Greenspan,  J.  M.  Ruohonicmi,  and 
M.  R.  Hairston,  The  interaction  of  a  magnetic  cloud  with  the  Earth:  Ionospheric  convection  in  the 
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northern  and  southern  hemispheres  for  a  wide  range  of  quasi— steady  interplanetary  magnetic  field 
conditions,  accepted  by  /  Geophys.  Res.  in  1992,  to  be  published  in  199.". 


Papers  submitted 

Emery,  B.  A,,  D.  J.  Knipp,  A,  D.  Richmond,  N.  U.  Crooker,  M.  R.  Hairston,  ct  a!..  Ionospheric  re¬ 
sponse  to  an  interplanetary  magnetic  cloud  —  1.  Overview  of  analysis  lor  1 98s  january  14  - 1 5,  sun  - 
mittedto/.  Geophys.  Res.  1992. 

Knipp,  D.  J.,  B.  A.  Emery,  A.  D.  Richmond,  N.  U.  Crooker,  R.  A.  Hcelis,  M.  R.  Hairston,  J.  A.  Cum¬ 
nock.  et  ai..  Ionospheric  response  to  an  interplanetary  magnetic  cloud  —  II.  A  case  study  of  14 
January  1988,  submitted  to/.  Geophys.  Res.  1992. 

Rich.  F.  J.  and  M.  R.  Hairston,  Large  scale  convection  patterns  observed  by  DMSP,  submitted  to/. 
Geophys.  Res.  1992. 

Papers  in  final  stages  awaiting  submission 

Hairston,  M.  R..  R.  A.  Hcelis,  and  F.  J.  Rich,  Comparison  of  simultaneous  potential  distribution  pat  - 
terns  in  the  northern  and  southern  polar  ionosphere,  to  be  submitted  to/.  Geophys.  Res.  1993. 
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